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ABSTRACT 

Aims. We applied a decomposition method to the energy dependent pulse profiles of the accreting binary pulsar A 0535+26, in order 
to identify the contribution of the two magnetic poles of the neutron star and to obtain constraints on the geometry of the system and 
on the beam pattern. 

Methods. We analyzed pulse profiles obtained from RXTE observations in the X-ray regime. Basic assumptions of the method are 
that the asymmetry observed in the pulse profiles is caused by non-antipodal magnetic poles and that the emission regions have ax- 
isymmetric beam patterns. 

Results. Constraints on the geometry of the pulsar and a possible solution of the beam pattern are given. We interpreted the recon- 
structed beam pattern in terms of a geometrical model of a hollow column plus a halo of scattered radiation on the neutron star surface, 
which includes relativistic light deflection. 
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1. Introduction 

The Be/ X-ray binary A 0535+2 6 was discovered by ArieJ V 
in 1975 (Rosenbe rg et al.l Tl975l) during a giant outburst. The 
system consists of a neutron star orbiting the optical com- 
panion HDE 245770 on an eccentric orbit (e = 0.47 ± 0.02 ) 
of orbital period P orb = 111.1 ± 0.3 d (Fin ger et "all l2006l) . 
The source presents quiescent X-ray emission with a luminos- 
ity of Lx < 10 35 36 ergs _1 , sometimes interrupted by "nor- 
mal" outbursts (Lx ~ 10 37 ergs -1 ) linked to the periastron pas- 
sages of the neutron star, and less frequent "giant" outbursts 
(Lx > 10 37 ergs~') of longer duration and less clearly re- 
lated to the orbital phase (see. e.g..lGiovannelli & Graziatill992t 
Kendziorra et aD ll994t iFinger et al.lll996l) . The system presents 
two cyclotron resonant scattering features at E ~ 45keV and 
E ~ 100 keV, from which a magnetic field strength of B ~ 
4 x 10 12 G is inferred (iKendziorra et aDl994|lGr"ove "et al.lll995t 
iKretschmar et al.ll2005HCaballero et alj|2007l) . 

Pulsations are observed to have a period of f sp i n ~103 s, typi- 
cally with spin-up during stronger outbursts and spin-down dur- 
ing quiescent periods^. The pulse profile evolves from a com- 
plex profile at lower energies to a simpler, two-peaked structure 
at higher energies. This behavior is observed in several accret- 
ing X -ray pulsars. Similar to other sources (e.g.. IStaubert et all 
1980), individual pulses show strong pulse-to-pulse variations, 
while the average pulse profile is rath er stable, with slower v ari- 
ations over the course of an outburst (Caball ero et al.l l2008b). 

The basic concept of pulsed emission is well understood. 
Pulsed emission originates in regions close to the magnetic poles 



of a rotating neutron star with the magnetic axis misaligned 
with respect to the rotation axis. In contrast, physical model- 
ing of the pulsed emission turns out to be a complex task. Many 
processes are in fact involved in modeling pulse profiles, from 
the modeling of the emission regions and their local emission 
pattern to the formation of the pulse profiles seen by a dis- 
tant observer. Comparison of model ca lculations with observa - 
tions has been performed for in stance bvlWang & Welter! d 198 lb . 
iMeszaros & Nagell (1 19851) . and lLeahvl (Il99ll) 

A proper model calculation should include relativistic light 
deflection, which has a significant effect on the pulse shap^l 
(Riffert & Meszarosl ll988l) . For slowly rotating neutron stars, 
the metric around a neutron st ar can be approximat ed by the 
Schwarzschild metric (see, e.g., iPechenick et al.lU983l) . Due to 
the strong gravitational field around the neutron star, the X-rays 
will be observed at red-shifted energies. Geometrical models of 
filled and hollow accretion columns of accreting neutron stars, 
includ ing relativistic light de flection, were computed in iKrausI 
(2001) and lKraus et al.l (12003). These models give the beam pat- 
tern or energy-dependent flux of one emission region as a func- 
tion of the angle, as seen by a distant observer. Introducing the 
rotation of the pulsar and its geometry, i.e., the orientation of the 
rotation axis with respect to the direction of observation and the 
location of the two poles, the pulsed emission from each of the 
two poles (single-pole pulse profiles) observed by a distant ob- 
server can be modeled. The sum of the single-pole contributions 
gives the total pulse profile. 



1 see, e.g., the results of Fermi-GBM monitoring at 
http : //gammaray . nsstc . nasa . gov/ gbm/science/pulsars/ 



The importance of relativistic light deflec- 
tion in model calculations can be visualized in 
http : //www. spacetimetravel . org/xpulsarQ6/xpulsar86. html 
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An alternative method of analyzing pulse profiles is to 
start from the observed pulse profiles and, based on symme- 
try considerations, decompose the pulse profile into single- 
pole contributions. This is then transformed into the visible 
section of the beam pattern. This method has been success- 
fully applied to the accretin g X-ray binary pulsars CenX-3 
HerX-1 and EXO 2030+ 375 feraus et al.lll996l iBlum & Krausl 
2000, ISasaki et alj 1201 Ol respectively) and is applied in this 
work to A 0535+26. Preliminary results were presented in 
ICaballero et all d2008ah . 

2. Method 

2.1. Assumptions 

A major goal of this work is to reconstruct the beam pattern of 
the neutron star in the binary system A 0535+26. By beam pat- 
tern we refer to the emission of one magnetic pole, as seen by a 
distant observer, as a function of the angle 6 between the direc- 
tion of observation and the magnetic axis (see Fig.[T). The basic 
assumption made in this work, which is often adopted in model 
calculations, is that the emission regions at the magnetic poles 
are axisymmetric. The magnetic axis is therefore a symmetry 
axis of the emission from the neutron star. Under this assump- 
tion, the beam pattern from one pole is only a function of angle 
0, and therefore the pulsed emission of each magnetic pole or 
single-pole pulse profile has to be symmetric. One of the sym- 
metry points will be the instant when the magnetic axis is closest 
to the line of sight, and the other symmetry point will be half a 
period later, with the magnetic axis pointing away from the ob- 
server. 

In the case of an ideal dipole field, the sum of the two 
single-pole pulse profiles will give a symmetric pulse profile be- 
cause the two single-pole pulse profiles have the same symmetry 
points. However, the observed pulse profiles of accreting X-ray 
pulsars are generally asymmetric. To explain the asymmetry ob- 
served in the pulse profiles, we assume a dipole magnetic field 
with a small offset from an ideal dipole field (see Fig.|2|. A small 
deviation of one of the magnetic poles from the antipodal posi- 
tion makes the symmetry points of the two single-pulse profiles 
different, causing the asymmetry in the total pulse profile. 

Another assumption made in this work is that the two emis- 
sion regions are the same, i.e., have the same beam pattern. This 
implies that each of the two poles will make one section of the 
same beam pattern visible. Depending on the geometry of the 
neutron star and the angle of observation, those two sections will 
have coincident parts in some cases. This assumption w as tested 
with the accreting pulsars Cen X-3 dKraus et alJll996l) and Her 
X-l (IBlum & Kr aus 2000). In these cases an overlapping region 
was found, in agreement with the assumption of two equal emis- 
sion regions. 

2.2. Decomposition into single-pole pulse profiles 

The first step of the analysis is to express the original pulse pro- 
file as a Fourier series. The total pulse profile F is written as 

F(<D) = -M + ["* cos(fc<D) + v*sin(M>)] + m„ /2 cos(%) (1) 

k=l 

where n is the number of bins of the original pulse profile, O the 
phase, and u k , Vk are given by 

i r* 

u k = - \ F(<D)cos(M>)«i<I> (2) 



i r +n 

v k = - F(d>) sin(/tO)£/<D (3) 
n J-n 

Equation[T]gives a valid representation of the original pulse pro- 
file at all phases if the Fourier transform of F approaches zero as 
the frequency approaches n/2. 

The single-pole pulse profiles f\ and fa are described by the 
following symmetric functions: 

n/2 

M<t>)=-c + J]c k co S [k(®-® l )] (4) 

k=i 

, n/2 

=~do + Y. dk cos{£[0> - (<D 2 + jr)]}, (5) 

L k=\ 

where Oi and <I>2 are the symmetry points of /i(O) and fa(Q>), 
respectively. Formally, a decomposition of F into two symmetric 
functions exists for every choice of their symmetry points <t>i and 
<I>2. For convenience, we use the parameter A := n - (<L>i - Q>2), 
which represents the azimuthal displacement of one pole with 
respect to the antipodal position (see Sec. 12.2. II ). All formal de- 
compositions will be contained in the parameter space <5i - A, 
with < Oi < n and < A < tc/2. Once the formal decomposi- 
tions are calculated, the following criteria are applied to reduce 
the number of decompositions to physically meaningful ones: 

- positive criterion: both /i(4>) and fai®) must be positive, 
since they represent photon fluxes; 

- non-ripple criterion : the single -pole contributions /i(<t>) and 
fa(<&) should not be much more complicated than the original 
pulse profile. Individual pulse profiles with many peaks that 
cancel out in the sum are not accepted; 

- the same symmetry points i>i and <I>2 must give valid decom- 
positions in all energy bands. 

Once a possible decomposition is found, the symmetry points for 
each of the two poles cDi and <t>2, and the parameter A, related 
to the position of the emission regions on the neutron star, are 
determined. 

2.2.1. From single-pole pulse profiles to beam pattern 

In Fig. Q]a schematic view of a rotating neutron star is shown. 
A spherical coordinate system is used with the rotation axis as 
polar axis. As explained above, the beam pattern is assumed to 
be axisymmetric, hence to only depend on the angle between the 
direction of observation and the magnetic axis 9. The value of 
changes with the rotation angle €>. Depending on the position of 
the poles with respect to the rotation axis and depending on the 
direction of observation with respect to the magnetic axis, we 
only observe a section of the beam pattern for each pole. 

Applying the cosine formula to the spherical triangle in 
Fig.Q] we obtain 6 as a function of the phase O: 

cos 8 = cos ©o cos 0, + sin ©o sin ©,■ cos(0 - <&,) (6) 

where ©o is the polar angle of the direction of observation, 0, 
the polar angle of the ith pole, and O,- one symmetry point for 
the ith pole. 

The intrinsic pulsar geometry is shown in Fig. [2] A complete 
description of the pulsar can be given in terms of the polar angles 
©i and ©2, and the difference in their azimuthal angles <J>i -O2 = 
7i - A. The angular distance 5 between the location of the second 
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Fig. 1 . Schematic view of a rotating neutron star from 
iKraus et all d 19951) . A spherical coordinate system is used, with 
the rotation axis as polar axis. 0, is the polar angle of the ith 
pole, ©o is the polar angle of the direction of observation. The 
angle 9 between the magnetic pole and the direction of observa- 
tion changes with the rotation angle O. 



Fig. 2. Intrinsic geometry of the pulsar. With the rotation axis as 
polar axis, the magnetic poles are located at polar angles ©i and 
©2- The angular distance 6 between the second magnetic pole 
and the point that is antipodal to the first magneti c pole gives 
the de viation from an ideal dipole field. Figure from lKraus et al.l 
(Il995h . 



magnetic pole and the point that is antipodal to the first magnetic 
pole can be used as a measure for the deviation from an ideal 
dipole field. From Fig.[2]comes 



cos S - - cos ©2 cos ©i + sin ©2 sin ©i cos A. 



Considering the beam pattern as a function of cos 9 and the 
single-pole pulse profiles as functions of cos(0 - O,), there is 
no distortion between the two functions because the relation be- 
tween cos 9 and cos(0 - O,) is linear (Eq.[6j. Therefore, once we 
have the single-pole pulse profiles, by plotting them as a function 
of cos(0 - Of) we obtain two different sections of the same beam 
pattern. These two sections are different because we see a differ- 
ent interval of 9 for each of the two magnetic poles. Depending 
on the geometry, these two ranges of 9 can overlap, so there will 
be an overlapping region in the beam pattern. In this region, at 
an instant O, the first pole will be seen at an angle 9. At another 
instant O, the second pole will be seen at the same angle 9. We 
can use Eq.|6]to express the relation between O and O in terms 
of the geometric parameters 

cot©o(cos© 2 -cos©0 sin 2 ~ 
cos(O-Oi) = — + , cos(0-0 2 )(8) 



sin @i 



sin @i 



that we write as 



cos(0 - <b\) = a + frcos(0 - 2 ), b>0 



(9) 



where the parameter a represents the shift between the two 
single-pole pulse profiles. 

The sections that the single-pole pulse profiles have in com- 
mon can be found by representing them as a function of cos(<E> - 
O,), so the values of a and b can be determined. Once this is 
done, both sections of the beam pattern can be represented as a 
function of the same variable q, defined as 



q := 



cos 9 - cos ©o cos 



sin ©o sin ©i 



(10) 



Using Eqs.|6]and|9) 

cos(0 - Oi) = q, 



(7) cos(6 - 2 ) = {q- a)/b. 



(11) 
(12) 



Since the relation between q and cos 9 is linear, by plotting 
the two single-pole pulse profiles as a function of q we obtain 
the total visible part of the beam pattern without distortion. If 
the values of a and b have been determined, by comparing Eqs.|8] 
and [9] we obtain two equations relating the three geometric pa- 
rameters @o, @i, and ©2. They can be solved for ©i and ©2 as a 
function of ©q: 



tan©i = 



—2a tan ©o 



(a tan© ) 2 + b 2 -T 

b tan @i 

tan ©2 = . 

a tan@o tan@i + 1 



(13) 
(14) 



To obtain the location of the magnetic poles ©i and ©2, an 
independent determination of the direction of observation ©o is 
necessary. If ©o is known, the position of the two magnetic poles 
can be obtained from Eqs. Qj] and [14] The pulsar geometry is 
then completely determined: the location of the poles ©i and 
©2, their displacement from the antipodal position A, and the 
direction of observation ©o. Once the geometric parameters of 
the system are known, the single-pole pulse profiles can be ex- 
pressed as a function of 9 using Eq.|6] The observable part of the 
beam pattern is in thi s way comp l etely reconstructed. Further 
details can be found in lKraus et al.l (119951) . 

2.3. Application to A 0535+26 

Our analysis is based on RXTE-HEXTE observations of 
A 0535+26 during its August/September 2005 normal out- 
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Fig. 3. Example of A 0535+26 original pulse profiles used in the 
decomposition analysis, obtained from RXTE-HEXTE observa- 
tions performed during its August/September 2005 normal out- 
burst (MJD 53614.66), in the energy ranges (18.27-30.90)keV, 
(30.90-44.53)keV, (44.53-59.07) keV, (59.07-99.85)keV from 
top to bottom. The fluxes have been normalized to unity. Two 
pulse cycles are shown for clarity. 



burst. Pointed observations performed between MJD 53614.66- 
53615.6 were used (observation IDs 91085-01-01-03, -01-01- 
05 , -0 1 -02-00, -0 1 -02-0 1 , -0 1 -02-02, -0 1 -02-03 , for a total expo- 
sure time of ~20ks). We extracted RXTE-HEXTE background- 
substracted light curves in the energy ranges 18. 27-30. 90keV, 
30.90-44.53 keV, 44.53-59.07 keV, 59.07-99.85 keV. After ap- 
plying barycentr ic and orbital corr ection to the light curves 
(ephemeris from Finge r et al.l 120061) . we folded them with 3 
bins, using the best pulse period from ICaballero et aD (|2008b). 
An example of the energy-dependent pulse profiles is given in 
Fig. [3] The pulse profiles selected for the analysis correspond 
to the main part of the outburst, as those profiles remain stable 
durin g the outburst and al s o during historical obse r vations (see, 
e.g.. [Frontera et all [19851: iKendziorra et al l [19941; iFinger et al l 
1996). The flux during the observations was F ( 3_ 5 0)kev ~ (1.66 — 
1.85) x 10~ 8 ergcirT 2 s _I , giving a luminosity of L(3-50)k e v ~ 
(0.79 - 0.885) X 1 37 erg s -1 assuming a distance of d — 2kpc 
([Steele etal.lll998l) . 



(Stee le et al.l 11998). The flux was determined from PCA and 
HEXTEspectra, that we modeled with a power law with a high- 
energy cutoff plus two Gaussian absorption lines at ~ 45 and 
~ lOOkeV. Details of the spectral and timing analysis of the ob- 
servations can be found in ICaballero et all (120071 12008b). 



2.3.1. Search for acceptable decompositions 

The energy-dependent pulse profiles are expressed as Fourier se- 
ries using an FFT algorithm. Twenty Fourier coefficients out of 
the 30 sampled points are enough to properly describe the orig- 
inal profiles. The highest frequen cy terms are no t considered in 
order to avoid aliasing (see e.g. iPress et a"flll992l) . These func- 
tions are then written as the sum of two symmetric functions 
/l(O) and fi{<&). The condition F(<b)=f 1 (<S)+f 2 (<&) allows us 
to obtain the coefficients q and dk (Eqs. @] |5) as a function 



of Uk and (Eqs. [2] |3). To search for physically meaning- 
ful decompositions, the <Di - A parameter space is divided in 
1 x 1 ° boxes. All the formal decompositions are represented in 
this plane. First we apply the non-negative criterion. This crite- 
rion is more restrictive when the minimum of the pulse profile 
is close to zero. Knowledge of the background level is partic- 
ularly important in this step. A small amount of negative flux 
is allowed t o account for the ~ 1 % uncertainty of the HEXTE 
background dRofhschild et al.|[l998l) . The non-negative criterion 
reduces the number of acceptable decompositions considerably, 
which are mainly driven by the higher energy pulse profiles 
(E=59.1-99.8keV) that have a pulsed fraction of -79%. The 
result of applying this criterion, which requires each decompo- 
sition to be valid at all energy ranges, can be seen in the Q>\ - A 
plane in Fig. |4] (left). The black regions indicate where positive 
decompositions have been found. We limit the search of physi- 
cally acceptable decompositions to those regions. 

The non-ripples criterion is then applied. The original pulse 
profiles in the considered energy ranges have two main peaks. 
We expect the single-pole contributions not to have many more 
peaks than the original one. The number of peaks of the single- 
pole pulse profiles are counted. A quality function is defined as 
the inverse of the total number of peaks. With this method we 
obtain the decompositions sorted according to their quality func- 
tion. Higher values of the quality function correspond to decom- 
positions that are simpler and not much more complicated than 
the original pulse profile. 

For each decomposition, there is a certain region in the <Di- 
A plane that contains other qualitatively similar decompositions. 
To handle the large amount of possible decompositions that have 
to be studied, similar ones are grouped together into types. A de- 
composition belongs to a certain type if its square deviation with 
respect to the chosen representative decomposition of that type 
is smaller than a certain ;f 2 = 10~ 3 . Grouping them into types 
considerably reduces the number of decompositions that are left 
to study, since we only consider one representative of each type. 
For instance, in the 18.27 - 30.90 keV range, 1863 decomposi- 
tions are grouped into 286 types. Figure [4] (right) shows where 
the highest ranked profile representatives were found, grouped 
into types. We examined all these decompositions in different 
energy ranges, dividing the parameter space in five regions A, B, 
C, D, and E, shown in Fig. [4] (right). 

As discussed in detail below, by combining different energy 
ranges and different observations, we find the best decomposi- 
tion of the original pulse profiles for the parameters of Region 
A. They are shown in Fig. [5] with the remaining unmodulated 
flux available to distribute between the two symmetric functions 
that cannot be determined from the decomposition. The minima 
of the two symmetric functions have been shifted to zero, so 
that the sum of the two symmetric functions plus the unmod- 
ulated flux reproduces the original pulse profile. As expected, 
the single-pole pulse profiles at all energy ranges are not much 
more complicated than the sum. The symmetry points obtained 
for this best decomposition, calculated as the average from the 
different energy ranges, are given in Table [1] We take 5° as the 
estimate for the uncertainty in the symmetry points, which is the 
approximate range in region A for which the decompositions are 
similar. 

The decompositions in Region C are discarded because they 
present a very strong anti-correlation in the main peaks which 
seems artificial, not expected from two independent emission 
regions. The single-pole pulse profiles are generated indepen- 
dently, so they should not have features that match exactly and 
cancel out in the sum. Decompositions in Regions B and D 
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Fig. 4. Left: Result of applying the non-negative criterion: decompositions outside the black area are discarded. Right: Highest 
ranked decompositions according to the non-ripples criterion grouped in five regions A, B, C, D, and E. The different symbols 
represent different energy ranges: + (18.3-30.9)keV, x (30.9^14.5)keV , A (44.5-59. l)keV, □ (59.1-99.8)keV. 
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Fig. 5. Left: original pulse profiles described with a Fourier se- 
ries, normalized to unity in four different energy ranges. Right: 
best decomposition of the original pulse profiles in two sym- 
metric functions, corresponding to region A. The vertical lines 
indicate the ranges for the symmetry points for both single-pole 
pulse profiles, <t>i (dashed lines) and <t>2 (dash-dotted lines). The 
dotted horizontal lines represent the unmodulated flux left avail- 
able to distribute between the two functions. 



Table 1. Symmetry points and azimuthal displacement of one 
pole with respect to the antipodal position for the best decompo- 
sition of the A 0535+26 pulse profiles. 



Pole 1 


Pole 2 


A = 7T — (®[ — ® 2 ) 


O, = 72° ±5° 
3>i +tt = 252° ±5° 


H> 2 = 285° ± 5° 
C> 2 + n = 105° ± 5° 


A = 33° ± 5° 



have also been discarded, because the single-pole pulse profiles 
present an anti-correlation in many small features that cancel 
out in the sum, which is also not expected from two indepen- 
dent emission regions. In Region E we find the "second best" 
decomposition for <Di = 82° ± 5° and A = 63° ± 2°. The 
single-pole pulse profiles are not much more complicated than 



Fig. 6. Left: original pulse profiles described with a Fourier se- 
ries, normalized to unity in four different energy ranges. Right: 
"second best" decomposition of the original pulse profiles in two 
symmetric functions, corresponding to region E. Dashed and 
dash-dotted lines have the same meaning as in Fig. [5] 



the sum, and present an energy evolution similar to that of the 
total pulse profile. However, we favor the solution in Region A 
compared to Region E because the single-pole pulse profiles in 
Region E present a strong anti-correlation in the main features, 
not expected from a physical point of view. Another argument 
for rejecting decompositions in Regions C, D, and E is that they 
all present higher values of A. Under the assumption of slightly 
displaced magnetic poles, lower values of A are more likely to 
be real. T his was the case in the analysis of the accreting pulsar s 
CenX-3 dKraus et al.ll 19961) andHerX-1 dBlum & Krauslh oOO). 
where the best decompositions were found for low values of A. 
A further argument against decompositions in Regions B, C, and 
D emerges in the reconstruction of the beam pattern from the 
single-pole contributions (see Sect. 12.3.2b . This argument does 
not emerge for Region E, and therefore the beam pattern ob- 
tained for Region E is also discussed below. 
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Table 2. Geometrical parameters of A 0535+26 from the decom- 
position analysis. See Fig.[2]for the definition of the angles. 



Pole 1 


Pole 2 


offset from ideal dipole 


0, ~ 50° 


2 » 130° 


6-25° 



2.3.2. From single-pole pulse profiles to geometry and beam 
pattern 

In the previous section d2.3.11 > we obtained a decomposition of 
the pulse profiles of A 0535+26 in two single-pole components. 
We therefore have the pulse profile for each magnetic pole for 
a given energy range as a function of the rotation angle As 
shown in Sect. l2.2.Tl using Eq|6l by representing the single-pole 
pulse profiles as a function of cos((D-cD,), we obtain two sections 
of the same beam pattern. We performed this transformation to 
the single-pole pulse profiles of A 0535+26. In contrast to the 
accreting pulsars HerX-1 and Cen X-3, an overlapping region is 
not observed for A 0535+26. 

However, under the assumption of two identical emission re- 
gions, the two sections can almost be connected to each other, 
with a small gap in between in the beam pattern that remains 
unobservable to us. This means that, owing to the geometry of 
the neutron star and to its rotation, we are observing two differ- 
ent sections of the same beam pattern. However, unlike the case 
of an overlapping region, we cannot determine the parameters a 
and b of Eq. [9] from a fit. But we can estimate their values. We 
combined the sections of the beam pattern in all energy ranges 
using different values of a, which represents the shift between 
the two single-pole pulse profiles. The best estimate for the pa- 
rameter a is a - 2 .2 + 0.1. In the case of antipodal poles, b — 1 
dKraus etalJll995h . Since we are assuming a small distortion, b 
should be close to 1, so we make the assumption b = 1. 

With these estimates for a and b plus the direction of obser- 
vation ©o, it is possible to obtain the location of the poles ©i and 
©2. The angular distance between the location of the second pole 
and the point that is antipodal to the first pole 6 can be estimated 
using Eq. [7] In Fig.|7]the constraints on the geometry of the pul- 
sar are shown. The values of ©i, ©2 and 6 are represented for 
all possible values of the direction of observation ©q. The mass 
transfer from the optical companion is expected to align the ro- 
tation axes of the binary system on a shorter timescale than the 
system's lifetime. Assuming that the rotation axis of the neu- 
tron star is perpendicular to the orbital plane, the inclination of 
the system is identical t o the direction of observation i = ©o. 
iGiovannelli et al.l (|2007) determined an inclination for the sys- 
tem of i — 37 +2 °. We can therefore obtain the location of the 
poles ©i and ©2. Table[2]gives the estimated values we find for 
the position of the magnetic poles and the offset. 

It is then possible to plot the reconstructed sections of the 
beam pattern as a function of 0. The two sections of the beam 
pattern are reconstructed for e (13°, 87°) and e (93°, 167°). 
Figure [8] (left panels) shows the reconstructed beam pattern or 
emission of one magnetic pole as a function of 0, with = 0° 
meaning the distant observer is looking down onto the magnetic 
pole, and = 180° meaning the observer looks at the pole from 
the antipodal position. Figure [8] (right) shows the same beam 
pattern in polar representation. 

For regions B, C, and D, it was not possible to connect the 
two single pole contributions, providing a further argument for 
discarding them. For the "second best" decomposition (Region 
E), it has been possible to connect the reconstructed sections of 
the beam pattern in a similar way as for Region A, also suggest- 
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Fig. 7. Constraints on the pulsar geometry. The vertical lines 
indicate the orbital inclination (solid line) and its error range 
(dashed lines). 



ing that the two emission regions are the same. We can connect 
them with the same values of a and b as in Region A. This im- 
plies the same values for the polar angles of the magnetic poles 
@i and ©2. The value of A is higher than in region A, so that the 
offset from an ideal dipole field is also higher, 6 ~ 48 °. In Fig. [8] 
left (right panels) the beam patterns for the best decompositions 
in region E are represented for comparison. They differ slightly 
from those in Region A, but the main features of the two func- 
tions are similar. The single-pole pulse profiles also had very 
similar shapes, just different symmetry points. 

2.3.3. Interpretation and modeling of the beam pattern 

A characteristic feature of the reconstructed beam pattern is a 
minimum observed in the flux between w 30° - 40° (Fig. [8] 
left). This feature is present at all energies, leading us to believe 
that it could be related to the geometry of the accretion. Filled 
column models give a beam pattern in which the flux decreases 
at low values of 0, corresponding to the instant when the observer 
looks along the accretion stream dKraus et al.l2003l) . Introducing 
a hollow column plus a halo created on the neutron star surface 
around the column walls from scattered radiation emitted from 
the walls would explain the increase in flux at low values of 
and the minimum as increases when the observer looks directly 
into the column. 

Another characteristic feature of the reconstructed beam pat- 
terns at all energies is a steep increase in flux at high values 
of (0 > 120°, meaning the observer looks at one pole of 
the neutron star from the antipodal position). This can be due 
to gravitational light bending, making the emission from one 
pole visible from all directions and brightest from the antipo- 
dal position. Such a steep increase in flux at high values of has 
been observed in mod el calculations, p roducing a maximum in 
the flux at 9 = 180° dKraus et al.ll2003l) . This feature has also 
been observed in the reconstructed beam pa tterns of other X-ray 
pulsars (HerX-1 a nd EXO 2030+375, see iBlum & Krausl feOOO. 
ISasaki et alJl20Tob . 

To obtain estimates on the size of the accretion column, we 
applied a phenomenological model of a hollow column to the 
case of A 0535+26. The model has been computed for values of 
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Fig. 8. Left: Reconstructed beam pattern of A 0535+26 or emission of one magnetic pole as a function 8, as seen by a distant 
observer. The beam pattern is shown in linear representation for Regions A (left panels) and E (right panels) for different energy 
ranges. Right: reconstructed beam pattern for Region A in polar representation, for the energy ranges 18. 27-30. 90keV (solid line) 
and 59. 07-99. 85keV (dashed line). = 0° means that the distant observer is looking down onto the magnetic pole, and 8 = 180° 
means that the observer looks at the pole from the antipodal position. 



6 e [0,40]°, to try to reproduce the minimum observed in the 
reconstructed beam patterns of A 0535+26 for low values of 6 
(Fig. [8] left). The applied model describes medium-luminosity 
pulsars for which the local luminosity per pole is on the order of 

L,^L^M.5.10^(-f-)(-^-^. (15) 
4.y/2 cr s s \0.1rad/\1.4 M B ) cr s 

Here, a is the half opening angle of the base of the accretion 
funnel, Lgd = 4-jrcm p GM/a-j the Eddington luminosity, M the 
mass of the neutron star, cr x the Thomson scattering cross sec- 
tion, and cr s the magnetic scattering cross section averaged over 
direction and frequency. In this case a radiative shock is expected 
to occur at a height that is a small fraction of the stellar ra- 
dius. The shock separates a region of freely falling plasma above 
from a settling region of nearly stagnant plasma below, with ra- 
diation mainly originating in the shock and the settling region 
and escaping from the si des of the subshock column or mound 
dBasko & Sun vaevll 19761) . 

The modeling is performed as in lKraus et al. I J2001 and in- 
cludes a small accretion column, the formation of a halo on the 
neutron star surface, an accretion funnel delimited by magnetic 
field lines, magnetic scattering in the accretion stream, and rel- 
ativistic light deflection. This is an extension of the purel y ge- 
ometric hollow column models presented in lKraus | (|200ll). A 
detaile d study of this model will be presented in iRraus et alj 
(I2010h . 

Beam patterns were computed in the range 8 < 40° for ac- 
cretion onto a neutron star with the canonical values of mass 
M = 1.4 M and radius r n = 10 km. According to the observa- 
tions of A 0535+26, we take an asymptotic luminosity per pole 
of L M = 8 x 10 36 ergs _1 , corresponding to a local luminosity 
per pole of L = L 00 /(l - r s /r n ) = 1.4 x 10 37 ergs~' that is close 
to the limiting luminosity in Eq. [15] for typical values of a as 
studied below. The radiation is emitted from the side of the col- 
umn below a radiative shock set at r t = 10.5 km, and the local 
emission is taken to be an isotropic blackbody at the effective 
temperature of the column wall. This simple emission pattern is 
based on more detailed studies that have shown that the bulk of 



the radiation is expected to escape from the side of the accre- 
tion funnel (Wang and Frank 1981) and that, in static columns, 
the local emission in many circumstances is close to isotropic 
(Meszaros andNagel 1985). 

The inner and outer half opening angles a\ and a Q of the hol- 
low accretion funnel (see Fig. [9j» were varied as listed in Table[3] 
The table also lists the effective temperature of the column wall 
and the density that the accreting material has at the base of the 
free fall column. These parameters follow from the assumed ge- 
ometr y and the prescribed luminosity as described in lKraus et al.1 
(l2003h . The effective temperature of the column walls, when as- 
suming an isotropic black body, is derived from <xL^ ff A co = L, 
where <x is the Stefan-Boltzmann constant, A co the emitting area, 
and L the local luminosity per pole. The free-fall density in the 
plasma rest frame at the base of the free fall column (at r = r n ) 
is taken to be homogeneous over the polar cap, and derives from 
poc 3 /3y(y - l)A cap = L where P = v/c a nd y — 1/ yl - ft 2 are 
evaluated at r — r n (see lKraus et al ] l2010l for more details). 

The results of the modeled beam patterns in the 8 e [0°, 40°] 
range for different values of the inner and outer half opening an- 
gles, for a photon energy of E — 7.6 keV for a distant observer 
(E = lOkeV at the neutron star surface) are shown in Fig. [10] 
As in Fig. [8] 8 — 0° means the distant observer looks down onto 
the magnetic pole. The minimum in the computed beam patterns 
is caused by the passage of the accretion stream through the line 
of sight. These modeled beam patterns reproduce the shape of 
the reconstructed beam pattern of A 0535+26 for low values of 
8 well (Fig. [8). In Fig. [10] (left) computed beam patterns for 
different values of the opening angle are shown. The position 
of the minimum moves to higher values of 8 when the open- 
ing angle of the funnel increases. In Fig. [TO] (right) computed 
beam patterns for different values of the column thickness are 
shown. The width of the minimum increases for thicker columns. 
By comparing the reconstructed beam patterns of A 0535+26 in 
Fig. [8] (left) for low value of 8 with the modeled beam patterns 
in Fig. [TO] we estimate the half-opening angle of the accretion 
stream in the case of A 0535+26 to be a Q = 0.2 rad or ~ 11.5° 
and the wall thickness to be a - ct{ = 0.06 rad or ~ 3.4°. We 
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Fig. 9. Geometrical model of hollow column (Kraus 2001). 

used these estimates for the half-opening angle and for the wall 
thickness to investigate the energy dependence of the computed 
beam patterns. We calculated beam patterns for photon ener- 
gies E = lOkeV, E = 20keV, E = 30keV, and E = 40keV 
at the neutron star surface (or E - 7.6 keV, E - 15.2 keV, 
E = 22.8 keV, and E = 30.4 keV for a distant observer). The 
computed beam patterns for different photon energies are shown 
in Fig. QT| The passage of the accretion stream through the line 
of sight produces a distinct minimum that only weakly depends 
on the photon energy as expected from a geometric signature. 



Fig. 10. Beam pattern models for a hollow column with a halo 
of scattered radiation on the neutron star surface. 9 has the same 
meaning as in Fig. [8] The minimum is produced when observer 
looks in the direction of the accretion stream. Left: fixed column 
thickness a - a\ = 0.06 rad and different opening angles a - 
0.2rad (solid line), a = 0.15rad (dashed line), a - 0.1 rad 
(dotted line). Right: fixed outer opening angle a\ =0.1 rad and 
different column thickness a Q - a\ = 0.02 rad (solid line), a - 
ot\ = 0.04 rad (dashed line) and a - a{ = 0.06 rad (dotted line). 



Table 3. Model parameters used in the computation of beam pat- 
terns for A 0535+26, for# e [0°,40°]. 



model 


ffj (rad) 


a (rad) 


fc B r eff (keV) 


po (10- 5 g/cm- J ) 


1 


0.08 


0.1 


4.1 


16 


2 


0.06 


0.1 


4.1 


9 


3 


0.04 


0.1 


4.1 


6.8 


4 


0.09 


0.15 


3.7 


4 


5 


0.14 


0.2 


3.5 


2.8 



3. Summary and conclusions 

In this work, a decomposition analysis was applied to the 
A 0535+26 energy-dependent pulse profiles. A dipole magnetic 
field is assumed with axisymmetric emission regions. The asym- 
metry in the total pulse profiles is explained with a small offset 
from one of the magnetic poles from the antipodal position. We 
find a physically acceptable decomposition of the pulse profiles 
that allows us to extract information on the geometry of the pul- 
sar. We obtain 0[ « 50° and © 2 ~ 130° for the position of the 
magnetic poles, and an offset of 6 « 25°. 

The visible section of the beam pattern was reconstructed. A 
characteristic feature of the reconstructed beam pattern at all en- 
ergies is a minimum observed in the flux between 9 w 30° - 40°, 
where 9 is the angle between the direction of observation and 
the magnetic axis. This was interpreted in terms of a simple ge- 
ometrical model that includes relativistic light deflection. The 
model includes a hollow column emitting isotropically black 
body radiation, plus a thermal halo created on the neutron star 
surface around the column from scattered radiation emitted from 
the column walls. Another characteristic feature of the recon- 
structed beam pattern is a steep increase in flux at high values of 




eo eo 



Fig. 11. Computed beam patterns for different photon energies 
for an outer opening angle a Q = 0.2 rad and a column thickness 
a Q - or; = 0.06 rad. 9 has the same meaning as in Fig. [8] 



9 (9 > 120°). This could come from gravitational light bending, 
which produces a similar feature in model calculations. 

We performed model calculations for different column thick- 
nesses and opening angles, and we found the best estimates of 
the half-opening angle and column thickness to be a a = 0.2 rad, 
a Q - a\ = 0.06 rad. We would like to stress, however, that 
this model is simplified, and we do not claim it is true in all 
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details, but it does reproduce the basic shape of the energy- 
dependent reconstructed beam pattern of A 0535+26 for values 
of 6 < 40 °well. Computation of beam patterns at different en- 
ergies has revealed a weak dependence of the minimum and its 
depth with the energy, suggesting that the minimum is mainly 
an effect of the geometry of the system, produced when the ob- 
server looks directly onto the accretion stream. This weak energy 
dependence on the minimum is also found in the reconstructed 
beam patterns of A 0535+26. 
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